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I give a short review of the current status of lattice QCD calculations of the spectrum of 
charmonium. 

I. INTRODUCTION 

In recent years the numerical simulation of lattice quantum chromodynamics (QCD) has devel- 
oped into an enormously successful method for obtaining precise values of several Standard Model 
parameters. It is the only truly ah initio method available for doing nonperturbative studies of 
QCD. It provides a unified framework for treating mesons and baryons involving charm and bottom 
quarks, whether in a heavy-light or heavy-heavy configuration. 

For charm physics the primary phenomenological objectives of lattice QCD are to help in the 
discovery and characterization of excited states, including exotics, to determine decay constants [H 
[2] and form factors, and to calculate electromagnetic transition rates. Theoretical objectives include 
providing guidance for effective field theory. 

Since the last Charm conference there have been no new comprehensive numerical lattice studies 
of exotics and excited quarkonium states. So I will review the current status only briefly. There 
has been exciting progress, however, in the high precision determination of masses of several of the 
lowest lying states involving heavy quarks. Progress here is, of course, prerequisite to an accurate 
determination and characterization of excited states. 

Lattice simulations have their limitations. While static properties, including masses and matrix 
elements are relatively easy to calculate, scattering processes, including charmonium production, 
inclusive processes, and multihadronic (more than two-body) decays are very difficult. 

In addition to the usual challenges of extrapolating to the continuum and physical light quark 
masses, heavy quark simulations face special difficulties. The lattice spacing a introduces a mo- 
mentum cut off of order 1 /a, typically less than 2 — 3 GeV with today's lattices. As the quark mass 
M approaches the cut off, discretization errors grow. Such mass-related errors can be controlled 
by a suitable formulation of the Dirac action for the heavy quark. Among them are (1) a lattice 
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FIG. 1: Charmonium excitation spectrum in MeV from Ref. 0. Left: J states. Right states 
including exotics. Lattice results (blue) are compared with PDG values (black) and various quark model 
results (purple and orange). Color online. 

version of nonrelativistic QCD[3], which converges rather slowly for charm, but well for bottom, (2) 
the Fermilab action[l] and its improvements O |6], good for both charm and bottom, and (3) the 
highly improved staggered quark action (HISQ) with discretization errors of order 0{al{aMcf), 
good for charm, but not so good for bottom with today's lattices. 

II. EXCITATION SPECTRUM AND EXOTICS 

In 2008 Dudek et a/.|7j published results of a comprehensive lattice QCD study of the exci- 
tation spectrum of charmonium. In order to access excited states they introduced a large set of 
interpolating operators and used a "variational" method to determine their masses. 

How well does this work? Figure [l] (left) gives an impression. Consider the 1 channel, where 
the ground state J/ip has a very clean signal. Six states were found, but only the first three shown 
were unambiguous enough that the authors were willing to associate them with known levels, 
namely the J/ip, ipC^S) and ■0(35'). The lattice levels were generally higher than experiment and 
the discrepancies grew from 12 to 82 MeV. They were also higher than quark model determinations, 
as shown. In the PC = H — channel this study also produced a couple spin-exotic states as shown 
in Fig. [TJ See Ref. [7J for more states. 

These results are pioneering, but they are deficient for a number of reasons. Sea quark effects 
were omitted, the calculation was done at only one lattice spacing, so an extrapolation to the con- 
tinuum is not possible, and open charm states were ignored. The Hadron Spectrum Collaboration 
is remedying these shortcomings [S]. 
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FIG. 2: Left: Charmonium levels from Ref. [ilO^ based on splittings from the spin-averaged IS' level with 
quark masses determined from the Ds and light pseudoscalar mesons and scale determined from upsilon 
splittings. Results for four lattice spacings from 0.18 fm (red) to 0.09 fm (blue) are shown. Black lines are 
PDG values. Right: complete low-lying heavy and light meson spectrum from Ref. showing five masses 
used to fix four quark masses and the lattice scale and three masses predicted in advance of experiment. 



To include sea quarks and allow extrapolation to the continuum and to physical light quark 
masses requires generating gauge field configurations (ensembles) with various lattice spacings and 
light quark mass combinations. The MILC collaboration makes publicly available a large archive 
of such ensembles with lattice spacing ranging from 0.15 fm to 0.045 fm and light quark mass ratios 
iTT-ud/n^s ranging down to 0.05, compared with the physical value of approximately 0.037[9]. 

MILC gauge configurations are being used by the FNAL/MILC and HPQCD collaborations to 
study the quarkonium spectrum. The FNAL/MILC campaign currently uses Fermilab quarks for 
both charm and bottom and HPQCD uses HISQ charm quarks and NRQCD bottom quarks. Most 
of the FNAL/MILC results presented here are from last year's study[IO], which, like the HPQCD 
study, used a limited set of interpolating operators and looked at only the low-lying states. A new 
FNAL/MILC campaign currently under way is aimed at excited as well as ground states using a 
large set of interpolating operators, higher statistics, smaller lattice spacings, and more accurate 
heavy quark mass tuning. I will give a preview of the improvement we hope to obtain. 



III. LOW-LYING LEVELS 



Broadly speaking, lattice simulations do quite well reproducing the low- lying hadronic states. 
Figure [2] gives an overview of the charmonium spectrum of the FNAL/MILC collaboration (left) jlOj 
and the full gold-plated meson spectrum of the HPQCD collaboration[llJ (right). To appreciate 
the precision of the calculation, we turn to a finer mass scale in the remainder of this minireview. 
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FIG. 3: Hyperfine splitting in charmonium extrapolated to zero lattice spacing. PDG values are in black. 
Left: FNAL/MILC^Oj. Units are rl « 1.58 Gev-\ Right: HPQCD[T1]. 
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FIG. 4: Recent HPQCD results for the mass (left)il4j and (right) [H]. 



Historically, the hyperfine splitting of the IS" state has proven to be a stringent test of lattice 
methodology. In Fig. [s] recent results of the FNAL/MILC collaboration (left) |10j are compared 
with still more recent results of the HPQCD collaboration (right) [12] . The FNAL/MILC calculation 
was based on Fermilab quarks. The extrapolated value is 117(11) MeV. In this result annihilation 
effects are not included. The HPQCD result is based on HISQ charm quarks. The pink band 
indicates the full error budget (±2 MeV), and includes a small increase due to annihilation, as 
suggested by perturbation theory. Determining level shifts due to annihilation is difficult. A 
recent nonperturbative (lattice) result estimates, instead, a decrease of 1 — 4 MeV from charm 
annihilation|13] . Thus, apparently, we have reached the point where annihilation effects cannot be 
ignored. 

Further examples of achievable precision in gold-plated quantities are shown in Fig. |4j The Dg 
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FIG. 5: Left: 25* — 15 splitting from Ref. [TU] showing possible complications from open charm. Right: 
Sample FNAL/MILC preliminary result for the IS hyperfinc splitting from the new analysis with much 
reduced errors. 

mass is obtained to 3 MeV accuracy from the mass splitting M{cs) — ^M{r]c)^^, and the Be mass 
is obtained to 10 MeV accuracy from the spUtting M{Bc) - M{bb)/2 - M(7/c)/2[l5]. 

Going beyond gold-plated quantities, the FNAL/MILC simulation results for the spin-averaged 
25" — 15" splitting, shown in Fig. [5] (left) disa grees with the experimental value pOj. Could this be a 
result of complications from the nearby open charm threshold? Those simulations do not currently 
treat two-body scattering states. A recent calculation by Bali and Ehmann with a small basis set 
suggests that they may be important for even states well below threshold, but further study is 
needed [16]. 



IV. CONCLUSIONS AND OUTLOOK 



We have seen dramatic improvement in our ability to reproduce the masses of the lowest lying 
charmonium states. This has been brought about by the use of improved quark actions, higher 
statistics, accurate tuning of the heavy quark masses, smaller lattice spacing, and the availability of 
gauge field ensembles that support an extrapolation to physical sea quark masses and zero lattice 
spacing. New analysis campaigns with multiple interpolating operators are underway both by the 
Hadron Spectrum Collaboration and the Fermilab Lattice/MILC collaborations. A preview of the 
hyperfine splitting from this new analysis is shown in Fig. [s] (right), a considerable improvement 
over Fig. [3] (left). These campaigns should provide more accurate information about excited and 
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exotic states. Treating two-body scattering states will continue to be a long-term challenge. 



J. Simone, these proceedings. 
Hee Chang Na, these proceedings. 
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